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RELAXATION FROM STEADY STATES FAR FROM EQUILIBRIUM AND THE 
PERSISTENCE OF ANOMALOUS SHOCK BEHAVIOR IN WEAKLY IONIZED GASES 

ROBERT RUBINSTEIN* AND AARON H. AUSLENDER* 


Abstract. The decay of anomalous effects on shock waves in weakly ionized gases following plasma 
generator extinction has been measured in the anticipation that the decay time must correlate well with 
the relaxation time of the mechanism responsible for the anomalous effects. When the relaxation times 
cannot be measured directly, they are inferred theoretically, usually assuming that the initial state is nearly 
in thermal equilibrium. In this paper, it is demonstrated that relaxation from any steady state far from 
equilibrium, including the state of a weakly ionized gas, can proceed much more slowly than arguments 
based on relaxation from near equilibrium states might suggest. This result justifies a more careful analysis 
of the relaxation times in weakly ionized gases and suggests that although the experimental measurements of 
relaxation times did not lead to an unambiguous conclusion, this approach to understanding the anomalous 
effects may warrant further investigation. 

Key words, weakly ionized gases, drag reduction 

Subject classification. Physical Sciences 

1. Introduction. Renewed attention has recently been given to experimental reports of anomalous 
behavior of shock waves in weakly ionized gases [1]. In an attempt [2] to isolate the mechanism responsible 
for these effects, shocks were measured in a decaying plasma following extinction of the plasma generator. 
It was argued that the measured relaxation time of the anomalous effects must correlate well with the the 
relaxation time of their unknown cause; this measurement might therefore suggest the mechanism responsible 
for the anomalous effects, or perhaps even identify it unambiguously. 

Relaxation times cannot always be measured directly; in such cases, they are inferred theoretically, 
usually by assuming relaxation from a near-equilibrium state. But since a two- temperature weakly ionized 
gas in maintained in a steady state by energizing the electrons against the energy sink provided by heat 
conduction [3], such a weakly ionized gas is really in a steady state far from equilibrium. By analyzing 
representative examples, the present work shows that systems which are initially far from equilibrium relax 
to equilibrium much more slowly than systems which are initially near equilibrium. 

Retarded relaxation is demonstrated analytically for two systems: the Kats-Kontorovich steady state 
[4], a steady state maintained by sources and sinks of particles and energy in a gas of particles interacting 
through a power-law potential, and a nonequilibrium state of a gas of light particles diffusing in a gas of 
heavy particles [5]. In particular, relaxation of the Kats-Kontorovich steady state following extinction of the 
sources and sinks is shown to be algebraic in time instead of exponential: the relaxation rate is computed 
theoretically to be l/t where t is time. 

The existence of non-exponential relaxation in systems far from equilibrium suggests that relaxation 
times in weakly ionized gases should not be computed based on relaxation from a near- equilibrium state. 
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while the first author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA 
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This result justifies a more careful analysis of the relaxation times in weakly ionized gases and suggests that 
although the experimental measurements of relaxation times did not lead to an unambiguous conclusion, 
this approach to understanding the anomalous effects may warrant further investigation. 

2. Steady solutions of the kinetic equations. The //-theorem of statistical mechanics [5] demon- 
strates that the Maxwell- Boltzmann distribution is the unique long-time limit of the distribution function of 
any isolated many particle system. If the system is near but not in thermal equilibrium, then the Chapman- 
Enskog expansion shows [6] that the Boltzmann equation 

(21) I = 

can be approximated by the BGK equation 

In Eq. (2.1), / = /(r,p,f) is the single-particle distribution function and ft(/,/) is the Boltzmann collision 
operator [5], In Eq. (2.2), f eq denotes the Maxwell-Boltzmann distribution function; this equation states 
that a near-eqmlibrium distribution function relaxes exponentially to a Maxwellian, where the relaxation 
time u> 1 is given 6 for a hard-sphere gas by 


(2.3) 
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In Eq. (2.3), p is the mass density of the gas, m is the mass of a gas particle, d is the particle diameter, and 
T is the temperature; units in which the Boltzmann constant k = 1 are assumed. 

Analogous exponential relaxation is directly evident in the Fokker-Planck equation 5 for a gas of light 
particles diffusing in a gas of heavy particles, 


d£ = 1 d(p 2 s) 

dt p 2 dp 

(2 - 4) —*{&+%}- 

In Eq. (2.4), m and p are the mass and momentum of a light particle and T is the temperature of the 
gas of heavy particles. In this problem, the Maxwellian distribution is recovered by setting the flux s - 0; 
relaxation to this state from a nearby state occurs over a time of order ca —1 with 


Although the //-theorem determines the long-time limit of the Boltzmann equation, it does not dictate 
the time required to reach this limit. We will consider initial states from which relaxation to equilibrium is 
particularly slow. Both the Boltzmann equation [4] and Fokker-Planck equation admit solutions representing 
non-Maxwelhan steady states far from equilibrium: these are solutions of the steady forms of Eqs. (2.1) and 
(2.4) which prove to be singular at zero or infinite momentum. Thus, the non- Maxwellian solutions of Eq. 
(2.1) are not integrable over all momenta: it will be shown that the integral representing the total number 
of particles diverges at infinite momentum. Although this divergence might seem to be of merely formal 
interest, it is actually another significant consequence of the //-theorem, which rules out any steady solution 
of Eq. (2.1) other than the Maxwellian. Similar singularities will be shown to exist for the non-Maxwelhan 
solution of the Fokker-Planck equation Eq. (2.4). 
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It follows that these distributions only exist in a finite region of momentum space. If they are extended to 
be zero outside this region, the result is a solution of the inhomogeneous steady Boltzmann or Fokker-Planck 
equation 

(2.6) Sl(f,f) = -S 


(2.7) 


1 d(p 2 s) _ 
p 2 dp 


where the source term S is a sum of delta functions. 

Since the source S in Eqs. (2.6) and (2.7) vanishes except over a finite number of spheres in momentum 
space, the far from equilibrium distributions satisfy homogeneous equations everywhere else. From this 
viewpoint, the slow relaxation of these states is immediately suggested by the observation that initially, 
df jdt = 0 except where the source is non-zero. Relaxation of the non-equilibrium distribution requires that 
the effects of removing the source diffuse over all of momentum space; this process can be very slow. 

3. The Kats-Kontorovich steady state. The collision integral is given in terms of the transition 
probability I/(pi,p',pi,p), where p^p' are post-collision and pi,p are pre-collision momenta, by 

n (/, f) = ^J dpidp'dPiW + E' — Ei — E) x 

(3.1) <5(Pi +p' - Pi -p)f / (p'i,P , -Pi>P){/7J ~ ffi} 


where /, f\ are pre-collision and f’-f[ are post-collision distribution functions, 


/ = 

fi = /(r,Pi,t) 
/' = /(r,p',t) 

(3.2) /{ = /( r,p'i,f) 


and E is the kinetic energy 

(3.3) 




The conservation laws of particle mechanics imply that for any distribution function / , 


(3.4) 


I dpn(/,/) = 0 

J pdp ft(/, /) = 0 

I E(p)dpSl(f,f) = 0. 


For non-equilibrium f, non-zero fluxes in momentum space 

Jm{p') = -f dpft(fj) 

Jp>p' 

(3.5) J e (p') = - f dpE(p)Q(f,f) 

Jp>p' 


3 



are possible, where J m is a particle flux, and J e is an energy flux. Thus, if 


S(po) — A 0 5(p - p 0 ) + Bo6'(p - p 0 ), 


then the corresponding fluxes are 


Jm=4n{ ° 

{ AoPq - 2B 0 p% if p > po- 


j e _ I 0 if P < Po 

2 m 1 AqPq - 45 0 Pq if p > po- 


It follows that the source distribution 


2 J , 1 


5 ~ 4Tr { P Z 6{p - po) + IZ 6 '( p - po ') ~ - TT^P-Pi)} 


satisfies the equations 


(3.10) 


J dp S = 0 
J pdp 5 = 0 
J E(p)dp 5 = 0. 


In view of Eqs. (2.6) and (3.4), Eq. (3.10) implies that the source term 5 does not alter the total mass, 
momentum, or energy of the system; a steady state with 5 as a source term is therefore possible. The fluxes 


(3.11) 

Similarly, 


0 if P < Po 

Jm = ^ J if Po < P ^ Pi 
if P > Pi 

Je = 0 . 


__jt i y ^ 

W { ^ rf(P “ P ° ) + ^3^(P“P°)-^(P-Pl)-i^(P-Pl)} 


satisfies Eq. (3.10) and yields fluxes 


JO if p < po 
Je — \ J 9 if p0 < P < Pi 
[ 0 if p > pi 
J m = 0 . 


An alternative formulation [7] is to consider the Boltzmann equation with external forcing 


F -fp =«(/’/)• 
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At any discontinuity of the distribution function /, a delta function singularity occurs. The singular source 
distribution S could therefore be replaced by an appropriately chosen external force; however, we will not 

pursue this formulation further. 

Kats et al. [4] show that if U is homogeneous of degree p, 

(3.15) [/(Ap'^Ap'^p^Ap) = A^pi.p'.pi.p) 
then the collision integral Eq. (3.1) vanishes 4 for the power-law distribution 

(3.16) /( P) = Cm^d^pt 12 | J | 1/2 p- 7/2 -" /2 

corresponding to a steady state with constant mass flux J = J m - In Eq. (3.16), p» is defined by the 
requirement that 


(3.17) 


ff(p'nP'iPi.P) = (yV u o(P i.p'.Pi.P) 


with Uo homogeneous of degree zero. They also demonstrate that the distribution Eq. (3.16) is local, in 
the sense that the collision integral converges when evaluated on this distribution taken over all momenta, 

provided 


(3.18) 

and the power law exponent is in the range 


-3 < p < —1 


(3.19) 


— < — 2 . 


It follows from Eq. (3.19) that the integrals representing the total number of particles and the total 


energy 


N = j dpf (p) 

E= J E{p)dpf{p) 


(3.20) " = J 

both diverge at large momentum; therefore, the steady state distribution function Eq. (3.16) cannot exist 
for all momenta. The distribution function must be cut off at some finite value p = pi- If a low momentum 
cutoff p = po is assumed as well, so that Eq. (3.16) applies only when Po < P < Pi, then this trucated 
distribution with fluxes of the form Eq. (3.11) is the solution of the inhomogeneous Boltzmann equation 
with a source of the form Eq. (3.9). 

4. Relaxation of the Kats-Kontorovich steady state. We will investigate the relaxation of the 
steady state Eq. (3.16) to thermal equilibrium following removal of the sources and sinks. It will be 
convenient to work with an infinite power-law distribution in which the number of particles is infinite. For 
the Kats-Kontorovich distribution cut off at momentum p u the limiting Maxwellian distribution satisfies the 

inequality 

(4.1) f° — ^ J>3/2 

where C is a constant. The mass density and temperature scale in terms of the limiting momentum pi as 


p N ~p j 
T ~ E ~ Pi 


-(m+i)/2 


■b*-3)/2 
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therefore 


(4-3) fo~p ( r U]/4 - 

It follows that in the limit Pl - oo, the limiting Maxwellian distribution vanishes pointwise, although its 

integral over all momenta becomes infinite. The infinite region of integration makes this combination of 
limits possible. 

The imtial relaxation rate of the Kats-Kontorovich distribution is 
<4 ' 4) " = 

where the hmctional derivative is evaluated at the Kats-Kontorovich distribution. The absolute value of 
e Hurt is written as J„ to indicate that this is its initial value; the particle flux will become momentum 
dependent and will decay with time as the distribution relaxes to equilibrium. 

(2.1) E o q fth!'forl USSeStS ^ S6ek a Similarity SOlUti ° n ° f thG time - dcpcndcnt Boltzmann equation Eq. 

( 4 - 5 ) /( P, t) = Cm^d.- 1 J 0 1/2 p- 7/2 -/ i/2 fl i(_ l_p(M+i)/ 2 p -M/2 tJ i/ 2 ^ 


Introduce the substitution 




in the collision integral. The result is the equation 

i = ,«/!««> x 

+ ( 9 /)4/(^+l) e , _ ( 9i) 4/( P +l) ei _ ( g jV(M+l) e J x 
+ (g')2/(M+l) 0 / _ ( 9l) 2/(^+l) 0l _ (g) 2/(^ + l) 0 x 

UH) 2/ ^o' lt (<z') 2/(/1+1) o', ( 9i )V(m+x) 0 1i(?) 2/(m+i) 0 ] x 

[(9 , i)" (7+/1 ) /(1+m) ^ ( <?' 1 ) ( 9')" (7+ ^ ) / ( i + ^ ) ^(?') 

( 4 - 7 ) -( 9 i)-(7+M)/(l+m) 0 ( 9i )( 9 )-( 7+#i , /(1+m) ^)] 

where locahty property has been assumed in order to integrate over all momenta. Locality of the relaxing 
distribution will be demonstrated later for a restricted set of exponents p. 

If we write Eq. (4.7) simply as 

(4 ‘ 8) ^( 9 ) *= f 9i, qMdMJ) - <f>(qi)4>(q)]dq , 1 dq , dq 1 


then power counting shows that 


K(A^i , A q , A qi , A q) — A a K(gJ , q f , ^ , g) 


where 


(4.10) 


a = -3. 


This homogeneity immediately implies that the integral equation Eq. (4.8) admits a solution 


(4.11) 


4>{q) = c q l . 
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The relaxing distribution function has therefore the asymptotic form 


(4.12) 


f = 


J ~ t 


and the flux decays according to the power law 

(4-13) 

The conclusion that the distribution function relaxes algebraically with time is in sharp contrast to the 

exponential relaxation of a near-equlibrium distribution. 

These calculations all require the convergence of the relevant collision integrals Therefore, both the 
original Kats-Kontorovich distribution and the relaxing distribution Eq. (4.12) must c oca^ er 
Bef. 4 show that both distributions are local provided -3 < /r < -2; for the distributions w,th -2 < p < 1, 

the relaxing distribution is nonlocal, and the calculation leading to Eq. (4.12) is incomplete. 

5 Constant flux solution of the Fokker-Planck equation. The Fokker-Planck equation for a gas 
of light particles diffusing in a gas of heavy particles is given by Eq. (2.4). The thermal equilibrium solution 
is found by setting the flux s = 0, leading to the Maxwell- Boltzmann distribution 


(5.1) 


f eq = aexp(- 


2 mT 


A steady state far from equilibrium is defined by Eq. (2.4) by making the flux , nomsero; the corresponding 
distribution function is 


(5.2) 


/(p)=^cxp(-- T 




dp cM £v )+Af ' 


where A is any constant. 

For large momentum, this distribution function has the asymptotic expansion 


(5.3) 


f(p) 


(mT) 3(mT) 2 
p 3 P 5 


Evaluating Eqs. (3.20) at large momentum, 


N ~ logp p — * oc 


'P 


■ 00 . 


(5.4) 

Also, the second par, o, Eq. (2.4) implies that the flux s is singular when p = 0 Therefore, like .hr i Kats- 
Kontorovich steady state, the far from equilibrium distribution described by q. ( . ) o y exis s i 

region po < P < Pi and must be maintained by suitable sources and sinks. , . . . 

If the sources and sinks are removed, the distribution relaxes to Maxwellian. The relaxation „me Is 
simply the time required for the effect of shutting down the source, a, Po.P, to diffuse over all o momen bum 
space. Since the equation is linear, this time can be estimated immediately as (p, - Po) /B ms«ad » / 

as predicted by Eq. (2.5) for relaxation of a ncar-.quilibrium dkdribution. Thus, the time toquued I for he 
nonequilibrium distribution to relax to equilibrium depends on the momentum range of the caution 
and is generally much larger than the time ,0 relax a near-equilibrium distribution. Note also that, as m 
the Kats-Kontorovich distribution, the large momentum tail contain, considerable energy m view of 

divergence in Eq. (5.4). 


7 



In comparing the results of Sects. 4 and 5 it is annarent ftof <-u 

explained as consequences of temperature inhomogeneities in the experiments [81 feh^Tfchi ^ ^ 

the reported anomalies do not reveal any unexpected or remarkahl 1 \ M ’ fr ° m thlS VieW P° lnt > 

:rr rr ^ 
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